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Abstract. We present further observations of molecular gas in head-on collisions of spiral galaxies, this time 
of the CO(J = 1 -> 0) and CO(J = 2 ->- 1) lines in the UGC 813 - UGC 816 system. UGC 813/6 are only 
the second known example of head-on spiral-spiral collisions, the first example being the UGC 12914/5 pair. 
Strong CO emission is present in the bridge between UGC 813 and 816, unassociated with stellar emission, 
just as in UGC 12914/5. The CO emission from the UGC 813/6 bridge, not counting the emission from the 
galaxies themselves, is at least that of the entire Milky Way. Collisions of gas-rich spirals are really collisions 
between the interstellar media (ISMs) of the galaxies. We show that collisions between molecular clouds bring 
H2 into the bridge region. Although the dense clouds are ionized by the collisions, they cool and recombine 
very quickly and become molecular again even before the galactic disks separate. Because the clouds acquire 
an intermediate velocity post-collision, they are left in the bridge between the separating galaxies. The star 
formation efficiency appears low in the molecular clouds in the bridges. We speculate that the pre-stellar cores in 
the molecular clouds may expand during the cloud collisions, thus retarding future star formation. Because the 
ISM-ISM collisions discussed here require a very small impact parameter, they are rare among field spirals. In clus- 
ters, however, these collisions should be an important means of ejecting enriched gas from the inner parts of spirals. 

Key words. Galaxies: spiral - Galaxies: evolution - Galaxies: ISM - Galaxies: interaction - Galaxies: individual 
(UGC 813, UGC 816) 
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C3 I 1. Introduction UGC 12914/5 bridge was found to be nearly 5 times that 

of the entire Milky Way. Due to the appearance of the syn- 
A new class of colliding galaxies has recently come to , u . , , , . n 1 1 ,1 

00 chrotron brightness contours and the magnetic held, the 

be recognized - galaxy encounters in which the neu- ljr <ri 1 om a /e a • 1 1 a „ T a , 

to b J UGC 12914/5 system is also known as the Taffy galax- 

tral interstellar media (ISMs) 01 the systems hit each . „ Iinr , io/i; • v. ai. a- 

v 1 ,, . les . I he UGC 813/6 pair shares these properties and thus 
other in collisions with very small impact parameters, 1 r ,, , r , a- , • T ai • 1 

.. .... , , . is a member ol the class ol taffy galaxies. In this work, we 

irrespective of the relative inclinations of the galactic ^ , , , , , TT /^/-i m o m 1 -i a • 1 a e i 

. , -i, , -i . , show that the UGC 813/6 bridge contains a lot 01 molecu- 

planes. I hese systems are characterized by a bridge with , . rl TI nrnon»/t m, . fTTnr , 010 

v J J b lar gas, again like UGC 12914/5. The centers o f UGC 813 

strong synchrotron emission and abundant atomic hydro- joic • 1 i, j- r^/i n\ u FTu Z 71 

j O, • ' — — — -| 1 . J and 816 were previously observed m CO(1-0) by Zhu ct al. 

gen IjCondon et al. 1993, hereafter C93). Two galaxy pairs J-, QQ( J, 



have clearly undergone su ch a collision - the UGC 12914/5 T . - , 

, nn „ , — TIL r . rir l , TTrin 010 m \<-< j — "1I it is commonly believed that molecular clouds are too 

(C93, Jarrctt et al. 1999) and UGC 813/6 (Condon ct al. „ . . , , ., , 

\nnnrl i TTTt r,nZ\ • /-( , a • -a- i , small and dense t o hit each other [i.e . they have a low 

2002, also VV 769) pairs. Contrary to initial expecta- „,.. r n „ „ , n '"j . , „ , , 

, . . . 1.1 rilling factor, e.g. Jog & Solomon 1992) or to be affected 

tions (C93), extremely strong CO em ission was detected , ... . . . — — . . . 

,, , .i ..„„ . /r , . , l r , . — i llno^nl by collisions with more diffuse atomic or ionized hydro- 
throughout the UGC 12914/5 bridge (Braine et al. 2003; , , T . ,f r 

tt: ~ — Tlln^nh i ai a i r i S en clouds. In addition to presenting the new results tor 

Gao ct al. 2003), showing that a large quantity of molecu- ? T ^,^, „ ,„ , , , . , 
, . , ii t r a ai nr\ • • r ai UGC 813/6, the goal of this work is to show that indeed 
lar gas is present as well. In fact, the CO emission horn the ' ' ° , 
collisions of molecular clouds are plausible and even in- 
tend offprint requests to: Jonathan Braine, braine@obs.u- evitable in collisions of gas-rich galaxies. As the surface 
bordeauxf.fr filling factor of the molecular gas in the spiral disks in- 
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creases, a head-on collision of the galaxies becomes radi- 
cally more efficient at drawing dense gas out of the parent 
disks and into the region between the two galaxies which 
separate after collision. The fate of the gas is unknown and 
probably quite dependent on whether the galaxies merge 
post-collision or not. Such collisions are, however, likely 
the only means of forcing much of the gas out of the inner 
regions of spirals because tidal forces primarily affect the 
less bound external regions. 

The stars are not affected by collisions of the ISMs of 
the galaxies and of course stars do not collide with each 
other. However, the optical appearance can be perturbed 
by the tidal forces which affect both stars and gas. The 
slower the relative velocities of the galaxies in the collision, 
the stronger the effect of the tidal forces will be. 

The UGC 813/6 bridge system is actually part of a 
triplet, the third member being CGCG 551-011 about 50 
kpc away and at the systemic veloc ity of the UGC 813/6 
pair, about 5200 km s" 1 . Following ICondon et alJ ((2002) 
we use the optical convention {v = cz) for all velocities. 
The atomic hydrogen (Hi) column densities are extremely 
high in the bridge region, reaching about 3 x 10 21 cm~ 2 . 
Two Hi peaks, Hie and Hiw, are also present in the outer 
parts of the common Hi envelope and are probably tidal 
features. The western Hi peak coincides with a blue object. 
The eastern Hi peak has no obvious optical counterpart 
bu t is close to the optic al tidal tail. 

ICondon et all (l2002ll suggests that UGC 813 and UGC 
816 are undergoing their first collision (some 40 to 50 
Myrs ago) and are now separating at about 400 - 500 
km s _1 . Contrary to the UGC 12914/5 system, the UGC 
813/6 disks are rotating in the same direction. CGCG 
551-011 is also a spiral galaxy but has been stripped of 
most of its Hi and has a strong compact radio continuum 
source. Without numerical simulations, we cannot deter- 
mine whether the tidal features could have been generated 
by a collision between UGC 816 and CGCG 551-011 some 
200 Myr ago, stripping CGCG 551-011 of most of its Hi 
and feeding the central source. 

After describing the observations of the UGC 813/6 
system, we present t he CO spectra and compare these 
with the Hi data from lCondon et all (|2002^l . Using the in- 
formation available for both the UGC 12914/5 and UGC 
813/6 pairs, we then show that molecular clouds can col- 
lide and drive large amounts of dense gas out of the spiral 
disks, despite earlier assumptions to the contrary. The dif- 
ferences with respect to previous work are explained and 
a general scenario proposed for how the gas, particularly 
the dense gas, gets into the bridge region for this type of 
collision. Taking Ho = 70 km s _1 Mpc -1 , we assume a 
distance to the UGC 813/6 group of 75 Mpc. 

2. Observations 

The UGC 813/6 bridge, galaxy centers, and the outlying 
Hi maxima were observed in the CO(l— >0) and CO(2— >1) 
to search for molecular gas. The observations were carried 
out with the 30 meter millimeter-wave telescope on Pico 



Veleta (Spain) run by the Institut de RadioAstronomie 
Millimetrique (PRAM) in August 2003. The CO(1->0) and 
CO(2— s-1) transitions at 115 and 230 GHz respectively 
were observed simultaneously and in both polarizations. 
A bandwidth of over 1300 km s _1 was available in both 
transitions using the two 512 x 1 MHz filterbanks at 3mm 
and the two 256 x 4 MHz filterbanks at 1mm. The band 
was centered around cz — 5200 km s _1 , where z is the red- 
shift, corresponding to redshifted frequencies of 113.3059 
and 226.6074 GHz. 

System temperatures were typically 200 - 300 K at 
3mm and twice as high for the CO (2-1) transition (T^ 
scale). The forward (main beam) efficiencies at Pico Veleta 
are currently estimated at 0.95 (0.74) at 115 GHz and 0.91 
(0.54) at 230 GHz. At the redshifted frequencies the half- 
power beamwidths are 22" and 11". All observations were 
done in wobbler-switching mode, usually with a throw of 
100" but sometimes more or less depending on the posi- 
tion observed, in order to be sure not to have emission in 
the reference beam. Pointing was checked on the nearby 
quasar 0133+476 every 60 - 80 minutes with pointing er- 
rors usually less than 3". 

Data reduction was straightforward. After eliminating 
the few obviously bad spectra or bad channels, the spec- 
tra for each position and each transition were summed. 
Only zero-order baselines {i.e. continuum levels) were sub- 
tracted to obtain the final spectra. Where no CO line was 
obvious, the line windows were based on the Hi spectra. 

3. Results 

The centers of both galaxies and 4 of the 5 bridge re- 
gions were clearly detected in both transitions. In Fig. 1 
we show the CO(1-0) and CO(2-l) spectra for the galaxy 
centers and bridge along with the Hi spectra at those po- 
sitions at 16" resolution, intermediate between the CO(l- 
0) and CO (2-1) angular resolutions. The northernmost 
bridge position shows positive flux at the appropriate ve- 
locities. No line is obvious in CO (2-1) although it covers 
an optically brighter region (a tidal arm) than the bridge 
center (Fig. 2). All of the pointing centers are shown in 
Fig. 2. In Fig. 3 we show the CO(1-0), CO(2-l), and Hi 
spectra of the eastern and western Hi peaks. No CO emis- 
sion was detected despite the strong Hi. 

Comparison of the CO (1-0) o bservations at 22" (this 
paper) and 55" (|Zhu et alJIliffll suggests that in UGC 
813 the molecular gas is more centrally concentrated than 
in UGC 816 because in UGC 813 the line is much brighter 
when observed at high resolution. 

In Table 1 we give the CO intensities for each transition 
and the corresponding H2 and Hi column densities for each 
of the positions observed. While the CO emission from the 
galaxies is stronger than in the bridge, the CO emission 
from the bridge is at least that of the Milky Way. For 
a CO - H 2 conversion factor of N(R 2 )/I C o = 2 x 10 2 
molecules cm" - 2 (K kms" 1 )- 1 Ce.g. lDickman et alJll98^ . 
the bridge is about half molecular and half atomic gas 
- some 2 x 10 9 M Q of each. This is based on the central 
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Fig. 1. CO(l^O), CO(2^1), and Hi spectra at the positions observed in the UGC 813/6 galaxies and bridge, using 
respectively black solid (white solid in the bridge), red dashed (yellow dashed in the bridge), and green dotted lines. 
The circle indicates the size of the CO(2-l) beam. UGC 813 is the edge-on spiral to the lower right and UGC 816 is 
to the left, less inclined and at a higher velocity. The asterisk marks a bright foreground star. For the CO spectra, 
the intensity scale is main beam antenna temperature in mK (-30 to 100 mK for the galaxy centers and -12 to 40 mK 
for the bridge positions). The Hi spectra are plotted such that equal CO(l— >0) and Hi intensity correspond to equal 
H-atom column densities in H2 and in Hi. It is thus apparent that in the bridge, the Hi and H2 column densities are 
similar, assuming the iV(H2)//co factor is correct. The underlying image is an R band image taken with the Instituto 
de Astroffsica de Andalucia (IAA) 1.5-m telescope on Pico Veleta. 



bridge position alone. For the other positions, it is difficult 
to separate the emission from the disks from that of the 
bridge, such that the above represents a lower limit to the 
true bridge CO emission. 

As can be readily seen from Fig. 1, the molecular and 
atomic gas line profiles are very similar. Furthermore, as- 
suming the iV(H2)/ico factor used is appropriate, the 
hydrogen column densities are also similar in the bridge 
region (the Hi spectra in Fig. 1 are scaled so that the 
same line area as in CO represents the same H-atom col- 
umn density) . The molecular gas dominates in the galactic 
centers and in the densest parts of the bridge. This is ob- 
viously not true of tidal tails, which preferentially bring 



material out of the Hi - dominated outer regions of spirals. 
The CO line intensities given in Table 1 for the bridge po- 
sitions are unheard of for tidal material. The two Hi peaks, 
at about 1-2 arcmin from the bridge and marked as Hie 
and Hiw, appear to be at the ends of tidal arms or tails 
and indeed the Hi emission is strong but no CO has been 
detected there so far (Fig. 3). The H2/H1 mass ratio in the 
tidal material is less than 10%, compatible with the less 
evolved Tidal Dwarf Galaxi es (TDG) for whic h CO de- 
tections have been obtained l)Braine et al.ll200l[) . Figure 4 
presents an Ha image of the UGC 813/6 system with Hi 
contours superposed. No star formation is observed near 
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Fig. 2. R ban d image taken with th e IAA 1.5meter telescope of the entire UGC 813/6 system with Hi column density- 
contours from ICondon et afl (J2002) . The centers of UGC 813 and 816 are marked by black triangles. The bridge 
positions observed are the five white triangles between the galaxies. These are the positions whose spectra are shown 
in Fig. 1. The triangles to the extreme West and East, slightly north of UGC 816, mark the positions of the Hi peaks 
whose spectra are shown in Fig. 3. Asterisks mark the positions of some galactic stars which could be confused with 
other features. The spiral galaxy seen nearly edge-on in the lower right is CGCG 551-011. 



the northeastern Hi peak (Hie), which is clearly not, or 
at least not yet, a TDG. 

To further distinguish between bridge and tidal fea- 
tures , we plot the synchrotron emission IjCondon et alJ 
|2002|) over the R band image in Fig. 5. Just as in the Taffy 
galaxies, the CO emission follows the synchrotron more 
than any tidal features. Strong extra-disk synchrotron 
emission is one of the characterizing features of ISM-ISM 
collisions and not of tidal arms or tails. 

In the Taffy galaxies (UGC 12914/5), the CO/Hi ratio 
in the bridge is higher. The CO and Hi line profiles in the 
UGC 12914/5 bridge are quite different; this is, however, 
probably a projection effect because we see the Hi from 
UGC 12914 as if it were in the b ridge while it is actually far 
behind (see iBraine et alJ2003l) . UGC 12914 and 12915 are 
generally bigger and more gas-rich than the UGC 813/6 



system. As we will see, this is important for the amount of 
molecular gas that can be stripped through cloud - cloud 
collisions. 

4. Why is H2 found in the bridges? 

In the following we will address the question of the ori- 
gin of the H2 found in the bridge between UGC 813 and 
UGC 816. As this galaxy pair is very similar to the Taffy 
galaxies UGC 12914/5 and we also found abundant molec- 
ular gas there, we will attempt to provide a more general 
explanation which applies to any Taffy-type interaction. 
We will assume that the 7V(H2)//co f actor adopt e d here 
is valid for the bridge molecular gas. iGao et all l|2003ft 
adopt a high er N(H?)/Ic.c> f actor in their work on UGC 
12914/5: IBraine et alJ l|20f)3h discuss why the factor they 
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Source 


RA 


Dec 


offset 


Ico(i-o) 


IcO(2-l) 


N(H a ) 


N(Hi) 




(J2000) 


(J2000) 


arcsec 


K km s" 1 


K km s" 1 


10 20 cm" 2 


10 20 cm 


UGC 813 


01:16:16.45 


46:44:25 


(0,0) 


12.2±0.9 


10.3±1.3 


24 


38 


UGC 816 


01:16:20.6 


46:44:52 


(0,0) 


11.1±0.6 


13.9 ±2.1 


22 


37 


Bridge 


01:16:18.4 


46:44:39 


(0,0) 


7.0±0.2 


3.7±0.4 


14 


35 


Bridge 






(7,-10) 


6.6±0.4 


7.4±0.5 


13 


27 


Bridge 






(-7,10) 


2.5±0.3 


0.7±0.4 


5 


32 


Bridge 






(10,7) 


6.7±0.4 


4.5±0.7 


13 


30 


Bridge 






(-10,-7) 


10.5±0.5 


8.3±0.6 


21 


35 


HlE 


01:16:24.31 


46:45:24 


(0,0) 


0±0.14 


0±0.4 


< 1 


22 


HlE 






(8,-8) 


0±0.10 


0±0.2 


< 0.6 


8 


Hiw 


01:16:06.1 


46:45:31 


(0,0) 


0±0.08 


0±0.2 


< 0.5 


10 



Table 1. Complete list of positions observed in CO, corresponding to the triangles in Figure 2. The offsets are in 
arseconds with respect to the (0,0) position of the source. Columns 3 and 4 give the CO(1-0) and CO(2-l) integrated 
intensities and cols. 5 and 6 the H2 and Hi column densities. The CO line window for the "Hie" positions is 5370 
to 5470 km s _1 and 5170 to 5270 km s _1 for the "Hiw" position, both based on the Hi spectra shown in Fig. 3. 
Uncertainties are la based on the rms noise and the line width. Upper limits to N(H2) are 3a. 




5000 5200 5400 5600 5000 5200 5400 5600 
Velocity (km/s) Velocity (km/s) 

Fig. 3. CO(1-0), CO(2-l) and Hi spectra of the "TDG 
candidate" positions. The brightness scale is for the CO(l- 
0) line and should be doubled for the CO(2-l) line (i.e. 
the CO(2-l) spectra have been divided by two). The Hi 
scale has been divided by 4 such that if the column den- 
sity of molecular gas were 25% of the atomic gas column 
density, then the spectra would be of the same height. It is 
immediately obvious that the CO intensities indicate low 
H2 column densities compared to Hi. The lower left panel 
shows the sum of the two eastern positions shown in the 
top two panels. The positions are shown in Fig. 2 where 
they are indicated by triangles. 



adopt is probably too high. The first question is perhaps 
whether the molecular gas could form from the atomic gas 
in the bridge. 



4.1. Could the H2 form out of the Hi in the bridge? 

The Hi column density is typically 2 — 4 x 10 21 cm~ 2 in 
the bridges. If the depth of the bridge is similar to its 
extent perpendicular to the collision direction, then the 
neutral gas is spread over some 10 - 20 kpc. The av erage 
volume density is then 0.1. iBraine et al.l l)200l[) esti- 
mate that 10 7 /n years (i.e. £ 10 8 yrs in the bridges) is 
necessary to convert 20% of the Hi into H2. Once some 
of the Hi has become H2 the process becomes less effi- 
cient. While local density enhancements could allow some 
H2 formation from the atomic gas, it is unlikely that the 
large quantities of H2 observed could have formed from Hi. 
H2 is in fact the dominant phase in the Taffy bridge and is 
roughly half of the UGC 813/6 bridge ISM. If the grains 
are warm then the ability of H-at oms to form H2 on the 
grain surfaces is greatly reduced l|Hollenbach fc Sakoeterl 
Il97ll) . Furthermore, the efficiency of Hi to H2 conversion 
is obviously reduced if grains are destroyed in the strong 
shocks generated in cloud collisions (whether in atomic or 
molecular clouds). 

There is a further reason for why clouds which are 
not initially dense will not be able to contract. In both 
systems, and particularly the Taffy galaxies, the spi- 
rals are mass i ve (se veral 10 11 M Q each from C93 and 
ICondon et alJ l|200 2)). As a result, the tidal forces they 
exert on clouds are s tronge r than in the Milky Way. 
Following iBania et alJ l|l986|) . for clouds to resist against 
tidal shear, their density must be 



200 ( 



3 kpc 
R 



) 2 nucleons cm 



where R is the distance from the cloud to the galactic 
centers. Being able to resist tidal forces is not sufficient 
for cloud survival but is certainly required. Thus, clouds 
which are not sufficiently dense after collision to form H2 
rapidly may not be dense enough to survive to form H2 
subsequently in the bridge. 
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4.2. GMC collisions as part of galaxy collisions 

The observations show that when spiral galaxies physi- 
cally collide, with the inner parts of their disks passing 
through each other, more than 10 9 Mq of molecular gas is 
found in the resulting bridge region between the two galax- 
ies. While tidal forces can drive gas out of spiral galaxies, 
they act on stars in the same way, such that it is diffi- 
cult to imagine how a bridge with so much gas but no 
stars could have been formed by tidal forces alone. I t has 
commonly been accepted (e.g. Ijog fc Solomonlll992() that 
Giant Molecular Clouds (GMCs) do not collide whereas 
the more diffuse clouds of atomic gas do. The goal of this 
section is to demonstrate that in gas-rich but not excep- 
tional spiral galaxies GMC-GMC collisions occur and are 
likely the source of the large quantities (1 - 10 xl0 9 M Q ) 
of molecular gas present in the bridges. 

In the calculations we will suppose t hat mo l ecular 
clouds are spherical and obey the so-called iLarsonl l)l98l|) 
relation with an average column density of N(H2) = 10 22 
cm" 2 . This is typical and corresponds exactly to what 
l.log fc Solomon! l|l992h assume. We define the following 
quantities: 

Ad, the projected cloud area, equal to 7ri? 2 j for identical 

spherical clouds where R c i is the cloud radius; 

n c i, the number density of such clouds; 

Met, the mass of a cloud of radius R c i; 

zh 2 , the half thickness of the layer of molecular clouds; 

Sh 2 > the- face-on H2 surface density of the spiral galaxies, 

expressed in Mq pc -2 , typically averaged over a telescope 

beam; 

£ c z, the H2 surface density of a cloud which is simply 160 
Mq pc" 2 for N(H 2 ) = 10 22 cm" 2 ; 

and the velocities V± , taken to be the encounter speed per- 
pendicular to the galactic planes and V|| , the relative cloud 
velocities in the planes (which are assumed parallel as in 
the Taffy galaxies). The link to extragalactic observables 
is provided by the integrated CO(l— >0) line intensity Ico, 
expressed in K km s , via the conversion factor taken to 
be N(H 2 )/I co = 2 x 10 20 H 2 cm" 2 per K km s" 1 . 

The mean free path A m f p of a spherical particle is 
A m fp = l/(4A c ;n c ;). For a gas of particles in random mo- 
tion, the mean free path is actually v2 shorter. The fac- 
tor 4 comes from the fact that clouds collide if they pass 
within 2R c i of each other. The number density of clouds 
is 

£h 2 .. 1 



2zh 2 



where M r i — S r ; x A r 



such that the mean free path is simply 



^mfp — 



2zh 2 S c 
4£h„ 



such that when Sjj 2 = S c ;, the mean free path of a cloud 
is 1/4 of the thickness of the molecular cloud layer 2zh 2 - 
This is in some ways a trivial result because £h 2 = £ c z 
means that the galactic planes are completely covered by 
a layer of molecular clouds so it seems obvious that they 
should collide when the disks pass through each other. 



CO observations show that many spirals have disk CO 
intensities of order Ico ~ 20 K km s _1 , corresponding to 
Eh 2 ~ 65 M Q pc -2 , such that the mean free path is less 
than the thickness of the molecular cloud layer. 

The mean free path formula above gives collisions even 
when only parts of the spherical clouds overlap. However, 
when there is a transverse velocity, such as in the Taffy 
galaxies where the counter-rotation yields high encounter 
speeds parallel to the galactic planes, then the path of a 
cloud in the molecular layer of the other system is in- 
creased from 2 zh 2 to 2 zh 2 



vl/v ± 



When the 



parallel and perpendicular velocities are roughly equal, 
this lengthens the path and increases the probability of 
cloud collisions by y/2. While the magnetic field energy in 
a cloud is greatly less than the kinetic energy of the colli- 
sion, magnetic fields are nonetheless quite strong in dense 
gas, with the field st rength increasin g as -^density, reach- 
ing about 100 fj, G l|Crutcherl lT999;l. As such, they could 
be expected to increase the effective cloud cross-section 
somewhat. If the galactic planes are not parallel and/or 
if the collision is not face-on then the cloud paths within 
the galactic planes are increased. 

The above discussion only considers collisions between 
molecular clouds, showing that rather than being rare, 
GMC - GMC collisions on large scales are inevitable in 
collisions of gas-rich spirals. After (inelastic) collision, the 
resulting cloud and possibly cloud fragments will have a 
velocity intermediate between the velocities of the galactic 
disks. They are thus left behind in the bridge region be- 
tween the separating galaxies. Whether or not the clouds 
fall back onto one of the galaxies depends on their velocity 
with respect to the galaxies and thus the collision velocity. 

4.3. Why previous results were different 

We choose to com pare here with the work of 
I.Tog fc Solomon! l(l992l) who explored this question and 
have been cited for why GMCs should not collide. 
Furthermore, the cloud properties they assume are in 
accord with our hypotheses - the only significant num- 
ber, the average H2 column density (or projected H2 
surface mass density) is the same as our adopted value. 
Ijog k, Solomon! l|l992j) define the mean free path by intro- 
ducing a volume filling factor for molecular clouds, which 
is the cloud volume multiplied by the number of clouds 
per unit volume in the molecular layer of a spiral galaxy. 
Hence 

/ = 4nR 3 cl n cl /3 

= (4rri&/3) x S H2 /(2zh 2 S c/ A c ,) 
= (4i? c ,/3) x £h 2 /(2zh 2 £ c 

after making the substitutions using the formulae above 
for identical spherical clouds (the hypothesis adopted by 
I.Tog fc Solomon! <|l992^ . Their mean free path is then 

A mfp = 2R cl /3f = 2R cl /3 x l/(47r/3)i^ n c ; = l/2A cl n cl 
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Fig. 4. Ha image after broadband continuum subtraction taken with the San Pedro Martir, Baja California (Mexico), 
2.12-m telescope with the same Hi column density contours as in Fig. 2. The dark vertical stripe through the bridge 
is due to the bright star above UGC 816. A fairly bright Hn region is seen in the southern part of the bridge. No Ha 
emission is observed near the NE Hi peak. 



which is twice the value we use. This may be to ensure that 
only collisions with comparable molecular column densi- 
ties are included, such that grazing encounters are entirely 
ex cluded. This is not the origin of the discrepancy. 

I.Tog fc Solomon! <)l992h specify the volume filling factor 
/ to be /gmc = 0-01 which we now express in terms closer 
to observables. Using the above definitions, 

E H2 = / x 2z H2 S c/ A c; x 3/(4vri?3 ; ) 

« J- x x ^M pc- 2 

0.01 130pc R cl 

« 6M Q pc" 2 

for 2z ho = 130pc, S. / = 16 0M w pc~ 2 , and R c i = 25pc, 
as m I.Tog fc Solomon! lll992l) . 6M©pc 2 corresponds to 
Ico = 2 K km s , which is weak for the inner parts 
of a spiral galaxy. This i s the origin of the discrepancy: 

calculations assume a very low 



molecular gas surface density. The value they assume is 
roughly the molecular gas surface density of the Milky 
Way averaged over the entire optical disk. In the molec- 
ular ring of our Galaxy the H 2 surface density is much 
higher. This is why grazing encounters should not be effi- 
cient at generating GMC-GMC collisions but collisions of 
the inner parts of spirals will be. 

5. The collision scenario 

Unfortunately, no simulations of colliding n eutral clouds, 
especi ally with magnetic fields, are available. lHarwit et alJ 
(1987) explore analytically what happens when two (big) 
disks of dense gas collide - we largely follow their scenario. 
The crossing or collision time of two clouds is 

r J~K~[ 500 km s- 1 
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Fig. 5. R band image with 1.4 GHz synchotron emission contours superposed. Triangles indicate the positions observed 
in CO and the asterisk a star. While the optical emission in the NW bridge position (near the northern tidal arm) is 
stronger than in the others, the CO emission is much weaker there. 



As two clouds collide, whether originally atomic or molec- 
ular, the colliding parts are instantly ionized and reach 
temperatures of T }t 10 6 K (depending on original colli- 
sion velocity). As the two clouds continue to move into 
each other, the central region will be compressed and will 
cool very rapidly. A similar scena rio is proposed in the 
case o f interstellar shocks in Fig. 1 of Hollcnb ach fc McKea 
l)l989|) . The pressure in the central region is very high so 
the gas will spr ead out perpendic ular to the collision direc- 
tion, as in the iMair et al.l l|l988j) simulations. How much 
gas flows to the side is highly dependent on the geometry 
of the clouds and of the collision. With the exception of 
gas flow towards the exterior, the c ollisions are concisel y 
desc ribed by Figu re 4 and Sect. 8 of lHarwit et all l|l987j) . 
The lStruckl l|l997(l simulations also suggest that ISM-ISM 
collisions can bring large amounts of gas into the bridge 
and that the gas is dense enough to cool, although he does 
not go into the details of cooling or the cloud collisions 
themselves. 



The hot gas (10 6 ~T < 10 7 K) cools very rapidly 
llCaetz fc Sa,lr ,eterl ll983t lHarwit et al]ll987h down to be- 
low 10 4 K. The cooling timescale, only a few years for 
molecular cloud collisions, is much shorter than the colli- 
sion timescale. Thus, the material at the barycenter of the 
collision will very quickly b ecome cool and neutral, pro- 
vided it is sufficiently dense. Hol lenbach fc McKed l|l989F) 
estimate the column density necessary to form H2 to be 
N1/2 ~ 4 x 10 21 cm~ 2 where the subscript 1/2 indicates 
changing half of the Hi into H2. The longest time scale 
is also the H2 formation time scale. All of the time scales 
vary as the inverse of the density so collisions of dense ma- 
terial have the best chances of producing predominantly 
cool molecular gas at the end of the collision. This means 
that the molecular gas we observe in the bridge region left 
the galactic disks in molecular form. 

Most collisions will obviously not be head-on collisions 
between identical clouds. Without detailed simulations of 
neutral cloud collisions, however, it is difficult to estimate, 
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for each geometry and possible cloud type (atomic, molec- 
ular), the fraction of the gas post-collision which will be in 
atomic, molecular, or ionized form. It is likely that head- 
on collisions of large atomic clouds will also produce (after 
ionization and cooling) a high fraction of H 2 . Given that 

(i) collisions between molecular clouds are not only pos- 
sible but inevitable on large scales once the surface filling 
factor (~ Eii 2 /Sci) of the colliding spirals is greater than 
about 0.2 (A m fp ~ thickness of molecular gas layer) and 

(ii) that collisions between dense clouds are most likely 
to generate cool dense (molecular) gas post-collision, it is 
reasonable to expect much of the H2 in the bridge regions 
to be due to collisions of molecular clouds. 

Observationally, if Hi - Hi cloud collisions were effi- 
cient in creating H 2 , then the Hi/H 2 ratio in the Taffy 
bridge would be relatively constant. However, the H 2 in 
the bridge decreases rather sharply from the center out- 
ward, just like the centrally concentrated CO emission in 
spirals, and the Hi is extended both in the bridge and in 
the spirals. The most likely explanation is that molecular 
cloud collisions generate much of the observed molecular 
gas in the bridges. 

6. Star Formation in the bridges? 

With such large quantities of molecular gas present, star 
formation rates of 1 - 10 Mq per year are expected in the 
bridges, if the gas is in a "normal" state, capable of form- 
ing stars as in spirals. From the Ha and CO emission of 
spirals, a star formation rate (SFR) of ab out lJlfe per year 
is estimated per 10 9 Mp> of molecular gas ijKennicutil lQQS: 
iBraine et al.ll200ll Sect. 3). The bridges are estimated to 
be about 20 - 30 Myr and 40 - 50 Myr ol d for the UGC 
12914/5 (T affy) and UGC 813/6 systems l|Condon et alJ 
h 9931120021) . Such times are well above the collapse times 
for cloud cores. 

From the CO emission from the bridges, the SFRs ex- 
pected in the UGC 12914 /5 and UGC 8 13/6 bridges are 



about 9 and 2 M Q yr _1 . barrett et a l. (1999) find that 



with the exception of the single giant Hn region near 
UGC 12915, the UGC 12914/5 bridge is devoid of stars 
and star formation. From t he slit spectrum of the UGC 
12914/5 bridge presented bv lBraine et all l)2003t) . there is 
no sign of star formation (except in the Hn region) al- 
thoug h the molecular gas surface densit y is nearly con- 
stant l|Braine et al.ll2003l iGao et al.ll2003h throughout the 
bridge. A single discrete Hn region is observed in the 
UGC 813/6 bridge (Fig. 4). The luminosity expected from 
bridges with several Moyr -1 of star formation is not ob- 
served. Because any stars formed in the bridge are nec- 
essarily quite young, could they still be dust-enshrouded, 

such that we would not see their light? 

From the Ha obs ervations ijBushousel Il987t 
iBushouse fc Wernerl Il990h of the Taffy pair, we esti- 
mate that the SFR is about 0.1 - 0.3 Mgyr -1 . From the 
Far-IR (IRAS) luminosity of the pair, about 4 x 10 Lq 
between 42 and 122//m, the SFR is between 5 and 10 
M yr _1 , more than half of which is in UGC 12915. In 



the Taf fy galaxies there is a lack of FIR emission in the 
bridge ijzink et alJ boOOj) but dust, perhaps only small 
grains, is detected via its Mid-IR emission. The flux in 
the bridge at 15/mi measured by l.larrett et ail l|l999h 
with ISOCAM is about 5% of that measured in each of 
the galaxies themselves. This is in rough agreement with 
the Ha based SFR so we estimate that the bridge SFR is 
between 0.1 and 0.4 M Q yr _1 . For the Taffy system, this 
is less than 10% of what is expected based on the CO 
luminosity. 

The cloud collisions and associated shocks are expected 
to destroy dust grains to a large extent. No observations 
relevant to this are available for the UGC 81 3/6 system. 
The s mall amount of dust in the Taffy bridge I Zink et all 
2000) means that it is unlikely that star formation is hid- 
den through extinction by dust. The dip in the FIR emis- 
sion in the bridge and the low SFR calculated from the 
MIR emission indicates that less star formation is occur- 
ring in the Taffy bridge than would be expected given th e 
strong CO emission. As discussed in lBraine et alJ (|2003j) . 
the A r (H 2 )/Jco ratio could be lower by a factor of a few. 
This would reduce the amount of "missing" star formation 
but even allowing for the uncertainty in the A(H 2 )//co 
conversion, the SFR still appears low. 

In the UGC 813/6 system, an Hll region is present in 
the bridge but contributes less than 2% of the Ha lumi- 
nosity of the system. The Ha emission surely suffers from 
extinction in the galaxies, particularly in UGC 813 which 
is seen edge-on. Although the pair has not been completely 
mapped in CO, we estimate that the bridge contributes 
between 12 and 30% of the global CO luminosity. Thus, 
again, the SFR is only some 10% of what would be ex- 
pected in "normal" molecular gas. Our tentative conclu- 
sion is that star formation in the bridge material has been 
retarded or even suppressed by the recent collision. 

In galaxies, stars form through the collapse of dense 
"pre-stellar" cores whose densities are above 10 5 cm -3 and 
column densities well above 10 22 cnr 2 . These cores are 
probably rather oblivious to what happens around them 
as they are gravitationally independent entities. Thus, if 
such cores survived into the bridge region, we would see 
the star formation occurring. 

To be ejected into the bridge region, however, the 
cores have experienced a sudden deceleration of several 
100 kms -1 . Just as for the rest of the cloud, the collision 
necessarily ionizes the dense cores as well. However, be- 
cause the probability of a dense core encountering another 
dense core is low, the dense core will continue deeper into 
the incoming cloud, ending up in a lower density medium 
than the barycenter of the two clouds. It may then be able 
to expand and lose its "pre-stellar" nature. The ionization 
and subsequent cooling of the clouds may result in some 
degree of homogenizing of the clouds, such that time is re- 
quired before overdensities can become pre-stellar cores. 

Is there any evidence for this in the line ratios? Only 
the Taffy system is strong enou gh in the CO lines t o have 
any 13 CO or HCN observations. IBraine et all l)2003h found 
that the 13 CO lines were quite weak in the Taffy bridge, 
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with line ratios of about 50 and 100 respectively in the (1- 
0) and (2-1) transitions respectively. Taken at face value, 
these line ratios indicate that the 12 CO lines are not as 
optically thick as they are in galactic clouds or in normal 
galaxies. In UGC 12915 the ratios are lower, with ratios of 
15 and 19 in the (1-0) and (2-1) transitions respectively, 
close to values in normal spirals. However, in IR-luminous 
galaxies with ve ry high SFRs. the 12 / 13 CO line ratios are 
also very high l|Casoli et al.l Il99^l and the star forma- 
tion is particularly efficient. In the IR-luminous galaxies, 
the HCN(l-O) l ine is strong, sometim es stronger than the 
13 CO(1-0) line llSolomon et all 19921) . In the Taffy bridge, 
we did not detect the HCN(l-O) line but the upper limit 
is only slightly below the 13 CO (1-0) line intensity in the 
bridge. Thus, the line intensities are compatible with a 
low opacity in the bridge region. 

7. Are head-on collisions cosmologically 
important? 

It is now generally accepted that galaxies were smaller at 
high rcdshift and that local universe galaxies (i.e. with 
small lookback times and hence current objects rather 
than early universe objects) are the product of collisions 
and accretion of smaller entities. The details of this pic- 
ture remain quite unclear, however, but these collisions are 
mostly tidal rather than the head-on collisions that are the 
subject of this work. We now consider two environments, 
field and dense clusters, to briefly explore whether ISM- 
ISM collisions could play a role in galaxy evolution for a 
large number of galaxies. 

To keep numbers simple, we will consider spiral galax- 
ies to have a radius of 10 kpc and take a dense cluster 
environment to have a density of 1000 galaxies per Mpc 3 
at a velocity of 1000 kms" 1 . The field galaxies are of the 
same size but have a density of 0.1 galaxies per Mpc 3 . 
Going through the same mean free path calculation as 
before, 

A mfp - l/(4A gal n gal ) = l/{4nR 2 n) « 800kpc 

but at 1000 kms -1 a galaxy moves 800 kpc in only 800 
Myr. In this extreme case, each dense cluster galaxy un- 
dergoes an ISM-ISM collision every Gyr! ISM-ISM col- 
lisions might then be a very efficient means of bringing 
gas out of galaxies and into the intracluster medium, es- 
pecially at earlier epochs when clusters were denser and 
galaxies more gaseous. The above numbers are probably 
overestimates: if the galaxies touch at a radius of lOkpc, 
the collision will not bring gas out of the inner parts - 
a more appropriate value of R would be 5 kpc; further- 
more, the galaxies should be assumed to have an average 
inclination of 45 degrees, reducing the area by V2; on 
the other hand, since the target galaxies move as well, 
the effective velocity is raised by about the same factor. 
Nonetheless, each cluster galaxy would still collide every 
few Gyrs, such that every cluster galaxy would undergo a 
head-on collision at least once and several times for some 



objects. Todays rich clusters are spiral-poor but they were 
certainly much more spiral-rich in the past, as well as be- 
ing denser. The gas ejected from the inner parts of galaxies 
would naturally be quite metal-enriched because star for- 
mation starts in the inner parts of spirals. 

A similar calculation shows that the situation is en- 
tirely different for field galaxies. At the much lower den- 
sity adopted for field spirals, only about one galaxy out 
of 1000 suffers this sort of collision in 10 Gyr. Even allow- 
ing for a density increase of n(z) w no(l + z) 3 , nothing 
changes out to a redshift of a few at least. We conclude 
that while these collisions are likely very important for 
galaxy clusters, they do not affect the evolution of field 
spirals. 

8. Conclusions 

Consistent with the Taffy galaxies, UGC 12914/5, the sec- 
ond such system, UGC 813/6, also contains large amounts 
of molecular gas in the bridge region between the galax- 
ies. UGC 813 and 816 are considerably smaller and less 
gas-rich than UGC 12914 and 12915, and the CO emis- 
sion from the UGC 813/6 bridge is about 4 times less than 
from the UGC 12914/5 bridge. Nonetheless, the CO emis- 
sion from the UGC 813/6 bridge is similar to that from 
the entire Milky Way and the H2 /Hi mass ratio is roughly 
unity. 

It is shown that for reasonably gas-rich, but not ex- 
traordinary, spiral galaxies a head-on collision will result 
in many collisions of molecular clouds. Some earlier calcu- 
lations had reached the opposite conclusion but they had 
supposed very gas-poor disks. Although the neutral (both 
atomic and molecular) clouds are ionized in collisions, the 
dense gas cools very rapidly allowing it to recombine and 
become molecular while still in the colliding disks. The 
new cloud is then left in the bridge because the collision 
has left it with an intermediate velocity. 

Some H2 may form through collisions of atomic gas 
clouds just as some of the originally molecular (or atomic) 
gas will remain ionized after collision. In addition to show- 
ing that GMC - GMC collisions occur in cases where the 
nuclei of the galaxies have passed within a few kpc of 
each other, several factors lead us to conclude that GMC 
- GMC collisions play a major role in bringing molecular 
gas into the bridg e regions. From a the o retica l point of 
view, the study bv iHollenbach fc McKed l|l989l) indicates 
that a column density of about 4 x 10 21 cm" 2 is neces- 
sary to form H2 . This is a large column density for purely 
atomic clouds. Furthermore, if Hi - Hi cloud collisions pro- 
duced the H2, the molecular gas would be as extended as 
the Hi in the bridge - this is not the case. Being much 
denser than atomic gas clouds, GMCs will be able to cool 
much more quickly after collision and ionization, forming 
the molecular gas we now observe in the bridge. 

The star formation in the bridge material is below that 
expected based on the molecular gas mass. We speculate 
that the cloud collisions may have eliminated many of the 
pre-stellar cores, such that in the few 10 7 years since the 
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collision the cores have not had time to reform and gen- Zhu, M., Seaquist, E. R., Davoust, E., Fraycr, D. T., & 
erate the strong stellar luminosity expected. Bushouse, H. A. 1999, AJ, 118, 145 

ISM-ISM collisions may provide an efficient mechanism Zink, E. C, Lester, D. F., Doppmann, G., & Harvey, P. M. 
for driving enriched gas out of the inner parts of spirals 2000, ApJS, 131, 413 
in galaxy clusters. Given that clusters were even denser in 
the past, most cluster galaxies have probably undergone 
such a collision at least once. Only a very small fraction 
of field spirals, even at high redshift, would have suffered 
a head-on collision. 

The future of the bridge material is as yet unclear. 
With a velocity of ^ 300 km s -1 with respect to one or 
another of the spirals, some of the gas is probably destined 
to fall back onto the parent spirals. Some however may 
remain in between the receding galaxies. 
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